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Abstract 

We describe an interpretation of convection in binary fluid mixtures as a superposition of thermal 
and solutal problems, with coupling due to advection and proportional to the separation parameter 
S. Many of the properties of binary fluid convection are then consequences of generic properties of 
2x2 matrices. The eigenvalues of 2 x 2 matrices varying continuously with a parameter r undergo 
either avoided crossing or complex coalescence, depending on the sign of the coupling (product of 
off-diagonal terms). We first consider the matrix governing the stability of the conductive state. 
When the thermal and solutal gradients act in concert {S > 0, avoided crossing), the growth rates 
of perturbations remain real and of either thermal or solutal type. In contrast, when the thermal 
and solutal gradients are of opposite signs {S < 0, complex coalescence), the growth rates become 
complex and are of mixed type. 

Surprisingly, the kinetic energy of nonlinear steady states is also governed by an eigenvalue prob- 
lem very similar to that governing the growth rates. More precisely, there is a quantitative analogy 
between the growth rates of the linear stability problem for infinite Prandtl number and the ampli- 
tudes of steady states of the minimal five-variable Veronis model for arbitrary Prandtl number. For 
positive S, avoided crossing leads to a distinction between low-amplitude solutal and high-amplitude 
thermal regimes. For negative S, the transition between real and complex eigenvalues leads to the 
creation of branches of finite amplitude, i.e. to saddle-node bifurcations. The codimension-two point 
at which the saddle-node bifurcations disappear, leading to a transition from subcritical to supercrit- 
ical pitchfork bifurcations, is exactly analogous to the Bogdanov codimension-two point at which the 
Hopf bifurcations disappear in the linear problem. 

PACS: 47.20.Ky, 47.20.-k, 47.20.Bp 

Keywords: binary fluids, double-diffusive convection 
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1 Introduction 



Convection due to two competing or cooperating effects is realized in a number of different physical 
systems: In the thermosolutal or thermohaline problem, vertical thermal and concentration gradients are 
both externally imposed. In convection in binary fluids with Soret effect, only the temperature gradient 
is imposed, but cross diffusion induces a concentration gradient with similar properties. An electrically 
conducting magnetic fluid may be subjected to a vertical or a horizontal magnetic field, the fluid layer 
may be rotated, or two solutes may be introduced. 

In the 1960s and 1970s, Veronis [0, 1 and other researchers |, |, |, |, 0, |, |, |l^, [l[ |l|l recognized the 
variety of behavior manifested by double-diffusive convection; comprehensive texts and reviews jl2 , 14 . Is] 



were written on the subject. One of the reasons for studying these double-diffusive systems is that all 
display a common basic set of phenomena. Both stationary and oscillatory instabilities occur, in other 
words pitchfork and Hopf bifurcations, the curves intersecting at parameter combinations which can 
be analytically calculated, at least approximately. In the 1980s, double-diffusive convection became the 
paradigm example in a renaissance in the study of bifurcation theory and dynamical systems, as attention 
was focused on precisely such intersections, re-interpreted as codimension-two points by Knobloch et al. 
^ and then by Brand et al. ||l|. 

There followed a divergence of efforts, roughly speaking between, on the one hand, detailed and rigor- 
ous mathematical analysis of the temporal complexity of thermosolutal convection and related problems 
in a small container, e.g. [0 |l8[ |l|, |48[ and, on the other hand, realistic and physical ex- 
ploration of the spatial complexity of the Soret and other problems in a large container for negative 

^ ^ HHHH H MMMMMMMMMM^ MMM^ S S 11 and 

positive 1^ ^,]2|,l3§ |9|, |g, gj, g |5^, |4| ^, |60|, |6l) values of the separation ratio. The first line 

of research has led to an understanding of the mechanisms producing global bifurcations, period-doubling, 
and chaos in two-dimensional convection in a confined geometry. The second line of research has resulted 
in bifurcation diagrams detailing the transitions in an extended periodic geometry, primarily between 
patterns of different spatio-temporal symmetry. For negative separation ratios, transitions occur between 
steady convection, standing waves, and traveling waves. For positive separation ratio, the transitions are 
between patterns of rolls and squares, and also between weakly and strongly convective regimes. 

In the present article, we will use an analytically tractable model of thermosolutal convection to 
examine the consequences of a simple idea: Consider a 2 x 2 real matrix 



(T-p 

7 



CTc 



Its eigenvalues are 




(1.1) 



(1.2) 



Each of the elements CTt, ctc and 7 depend on parameters r (scaled Rayleigh number) and S (separation 
parameter). We assume that, for each S, there is a value rint of r at which a-r and ctc coincide. We are 
interested in the behavior of the eigenvalues (t± in this vicinity. There are three possibilities, depending 
on the sign of P"f at r-mf If < 0, then cr_|^ and (t_ coalesce into a complex conjugate pair in the 
vicinity of rjnt- If f3j = 0, then a+ and cr_ intersect transversely. If f3j > 0, then (t_|_ and cr_ are real and 
cr+ > (J-, with a change of slope at r-mf This is the phenomenon known as avoided crossing. These three 
possibilities, shown in figure |l], correspond to the three scenarios observed when < 0, S" = 0, and S > 0, 
respectively. Moreover, we claim that this interpretation applies to the nonlinear steady states as well 
as to the linear stability problem. The remainder of the article is devoted to making these statements 
precise. 

The paper is organized as follows. Section 2 sets out the idealized free-slip thermosolutal problem that 
we will study. Although our calculations are limited to this analytically tractable case, we will present 
evidence throughout the paper that our conclusions may apply to double-diffusive convection problems 
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Figure 1: Eigenvalues of thermosolutal convection as a function of scaled Rayleigh number r for different 
signs of separation parameter S. S = —0.1: formation of complex conjugate pair. S ~ 0: eigenvalues 
cross transversely. S = 0.1: avoided crossing. 



in general. Section 3 gives the classic linear stability analysis and the bifurcations undergone by the 
system. The behavior of the eigenvalues is interpreted in the framework of avoided crossing vs. complex 
coalescence and also as primarily thermal vs. solutal. Section 4 analyzes the minimal five-mode Veronis 
model. We show that the steady states of the Veronis model obey a two-variable eigenvalue problem 
entirely analogous to the linear stability problem, with the energy playing the role of the eigenvalue. We 
present the consequences of this analogy, particularly for the saddle-node bifurcations and codimension- 
two points that occur for negative S. For positive S, the transition between Soret and Rayleigh regimes 
is analyzed in detail as exemplifying the phenomenon of avoided crossing. Section 5 sets out the simplest 
time-dependent system reproducing the primary bifurcations and steady states and, additionally, standing 
waves terminating in a heteroclinic orbit. However, the traveling waves that have been the subject of 
so much interest in double-diffusive convection are not accessible by our simplified approach; nor is the 
competition between patterns of rolls and squares. Section 6 is a brief summary and discussion. 



2 Thermosolutal problem 

We study the thermosolutal problem, in which thermal and concentration gradients are imposed, instead 
of the more widely studied and experimentally accessible Soret-driven problem, in which the concentra- 
tion gradient results from the temperature gradient, because its formulation is slightly simpler. The basic 
phenomena which we will investigate occur in both problems, which have in fact been shown by Knobloch 
[ p^ to be formally equivalent for the idealized boundary conditions we consider here: a two-dimensional 
geometry with free-slip upper and lower boundaries and imposed horizontal periodicity. Although all 
our calculations concern the idealized thermosolutal problem, wc will also refer to results in the liter- 
ature derived using rigid boundary conditions or concerning the Soret problem with various boundary 
conditions. Wc will not systematically refer to the related problems of rotating or magnctoconvcction 
(although these are also discussed i n Jl3| , |l^ , [T^ , [T^ , 48 ) nor to the equally vast experimental literature 



(for references, see, e.g. [50, p2|, p3|, p8| 



In the thermosolutal problem, the flow is driven by vertical thermal and concentration differences AT 
and AC imposed across a layer of height h. In the Boussinesq approximation, the density is assumed con- 
stant except in the buoyancy term, where it is taken to vary linearly with temperature and concentration, 
with thermal and solutal expansion coefficients pT and pc- There exists a steady motionless conductive 
solution, with temperature and concentration profiles depending linearly on the vertical coordinate z. 
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The stability of the conductive solution and the subsequent evolution of the system is determined 
by four nondimensional parameters. The three diffusivities - the thermal diffusivity k^, the momentum 
diffusivity or kinematic viscosity i/, and the solute diffusivity kc - are described by two nondimensional 
ratios, chosen to be the Prandtl number P =^ v/kt and the Lewis number L = kc/kt- For the linear 
problem, the Prandtl number is taken to be large; P = 0(1) or 0(10) suffices. The Lewis number is 
varied, but is considered to be small and, for the figures, will be fixed at the frequently-studied value 
of L = 0.01. The thermal Rayleigh number R = gpT^Th^ / {vkt) measures the imposed temperature 
gradient. The concentration gradient can be specified by an analogous solutal Rayleigh number, or by 
the separation parameter 5" = {pc ^C) / {pT ^T) , which may be either positive or negative. Length, 
temperature, concentration, and time are nondimensionalizcd by h, AT, AC, and h"^ / kt- 

It will be convenient to write the equations governing the evolution of the solution in terms of both 
the vertical velocity w and a velocity streamfunction cf), and temperature and concentration deviations T 
and C from their conductive profiles: 

dtf = V^f + w + ey • (V0 X Vf) (2.1a) 

dtC = LV^C + i& + Gy • (V0 X VC) (2.1b) 

dtV'^w = PV''w + Pi?a^(f + S'C)+ey9^(V(^x VV2(^) (2.1c) 

w = d^^ (2. Id) 

The boundary conditions are: 

f ^ip^dlip^ib^dlw at z = 0, 1 (2.2a) 

f , C, 4>, w 27r/A;-pcriodic in x (2.2b) 

where k may be fixed arbitrarily, or at the well-known value — tt/V2 which minimizes the convection 
threshold in the case of free-slip boundaries. 

Following, e.g. [|[ 0, we additionally scale time and, consequently, velocity, by: 

q"^ = k^ + TT^ (2.3a) 
and introduce the geometrically scaled Rayleigh number: 

r = Rk^/q^ (2.3b) 

Equations (2.1) become: 

dtf = q-^V^f + w + ey ■ {V^x Vf) (2.4a) 

dtC = Lq-^V^C + iv + ey{V^xVC) (2.4b) 

dtq-^V^w = Pq-'^V^w + Prk-^dl{f + SC) + q-^eyd^{y4>xVV^4)) (2.4c) 

w = dj (2.4d) 



3 Linear analysis 

3.1 Linear stability problem 

We begin by discussing the linear stability of the conductive solution, governed by: 

dtf = q-^V^f + w (3.1a) 

dtC = Lq-^W^C + w (3.1b) 

dtq-^V^w = Pq-'^W'^w + Prk-^dl{f + SC) (3.1c) 
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Solutions to (3.1) with boundary conditions (|2.2|) are of the form 



r(x, z, t) 

C{x,z,t) 
w{x, z, t) 



m 

C{t) I cos(fca;) sin(7rz) 
wit) 



(3.2a) 
(3.2b) 
(3.2c) 



The time dependence of the hnear system (p.l[) is 



T{t) \ [ T 

C(t) = exp(At) C 
w{t) j \ w 

wher e T, C , if a re sc alars denoting the amphtudes of the corresponding fields. Substituting (| 
and ( 2.3a ) into (3J_) yields the eigenvalue problem: 



(3.3) 

2D-(ra 



-1 














-i 


-I) 




(3.4) 


Pr 


P5r 









We simplify further by assuming that the Prandtl number P is infinite. Most quantitative results of 
interest to us depend only very weakly on P as long as P > 1, e.g. [Q. (Another approach to reducing 
the 3x3 problem to a 2 x 2 problem, which does not rely on P large, is given in Appendix 0.) The 
velocity amplitude w is then related to T and C by the algebraic equation: 



w = r(T + SC) 



(3.5) 



After eliminating the velocity via (|3.5D , the eigenvalue problem (3^) becomes: 

T \ ( r~\ rS \ f T 



A 



C 



rS-L 



C 



The matrix: 



r rS - L ) ~ \ Q -L 
M = Mo + rMi 



1 S 
1 S 



(3.6) 



(3.7) 



is of the form discussed in the introduction. The matrix Mq describes diffusion and the matrix Mi 



describes advection via (3.5). The diagonal elements of M: 



CTt = r — 1 
CTc = Sr ~ L 



(3.8a) 
(3.8b) 



can be viewed as the eigenvalues of a pure thermal and a pure solutal problem, with a coupling of: 

/37 = Sr^ (3.8c) 



The pure thermal convection problem for infinite Prandtl number is obtained from (3.6) by setting 
5 = 0. Bifurcation to steady thermal convection occurs when CTt = 0, i.e. at: 



r = = 1 



(3.9a) 



The pure solutal problem corresponds to convection driven exclusively by concentration gradients, 
i.e. incomplete mixing of the two species in the binary fluid. In the thermosolutal problem, the imposed 
concentration gradient is set to the value Sr, while in the Soret problem, the Soret coefficient S' is a 
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property of the fluid. Although we are studying the thermosolutal problem, it will be convenient for 
us to consider S as fixed and to vary the single control parameter r for all three problems: thermal, 
solutal, and thermosolutal. (In particular, all three problems have the same critical wavenumber for the 
idealized boundary conditions.) We therefore interpret ( 3.8b| ) as meaning that S is fixed and the onset 
of convection occurs at: 

r = rc = L/S (3.9b) 

There are two cases, depending on the sign of S (L, the ratio of two diffusion coefficients, is always 
positive): If 5 > 0, then adr) has positive slope and the conductive state is unstable for r > L/S > 0. 
If S* < 0, then <Jc{r) has negative slope and the conductive state is unstable for r < L/S < 0. 



3.2 Bifurcations 



We now return to the coupled thermosolutal system (^^). We begin by giving some exact results 
concerning the eigenvalues and bifurcations of (3^). Most of these results are well known, but we 
derive them here to illustrate our geometric interpretation and to prepare the analogy with the nonlinear 
problem of section ^. The results are summarized in the two large figures |^ and ^. 



The uncoupled thermal and solutal eigenvalues ( 3.8a ) and ( 3.8b ) intersect at 



1 - L 
1-S 



(3.10) 



The behavior of the thermosolutal eigenvalues of (3.6) near rint is determined by the sign of the coupling 
strength (product of the off-diagonal terms) S'r^, which is in turn determined by the sign of S. For 
S positive, the eigenvalues remain real for all r, with avoided crossing near rinf For S negative, the 
eigenvalues form a complex conjugate pair over an interval surrounding rint- 



The eigenvalues of (3.6) are: 



± 



(Tt — (Jc 



Sr^ 



= -[(l + ^)r-(l + L)] 

±^\/(l + S)^r^ - 2(1 - S)il ~ L)r + (1 - 
= /(r)±v/^ 



(3.11a) 

(3.11b) 
(3.11c) 



The linear function /(r) and the quadratic function g{r) are half the trace and a quarter of the discrim- 
inant, respectively, of the matrix M of (3.7). The qualitative behavior of (3.11) depends on whether g is 
positive for all r (S* > 0) or of both signs (5 < 0), and on whether g is a quadratic function of r {S ^ —1) 
or linear {S ~ —1). We write A± = cr± ± iuj to denote the real and imaginary parts of A, with (7+ = (t_ 
ff a; ^ 0. 



Appendix ^ derives properties of ( |3.1l| ) based on the matrices Mq and Mi in (3/7). Assuming that 
i S* 7^ — 1, then the discriminant Disci = (1 + S)"^ of matrix Mi is positive, and so the real parts cF±{r) of 
( |3.11 ) describe a hyperbola given by: 



S + L 
l + S 



{l + S)[r 



{l~L){l-S) 



S + L 
l + S 



S 



{i + sy 



and the imaginary parts +uj describe an ellipse given by: 



[l + Sf 



{l-L){l~S) 

{1 + sy 



-s 



{i + sY 



(3.12) 



(3.13) 



6 



(In the exceptional case S = —1, when /(r) is constant and g{r) is a linear function, then the curves 
a±(r) and ztu!{r) are parabolas; see figure 0.) 



The crucial quantity: 



A = S 



(1-L) 



(3.14) 



on the right-hand-side of both ( 3.12 ) and ( ^.13| ), an invariant for second-degree equations (see Appendix 
^ and ||6^ ) , distinguishes between avoided crossing and complex coalescence. The sign of A is determined 
by that of S. If S' > 0, the eigenvalues remain real: the branches (t_|_ and a- remain distinct and 
continuous over all values of r. If < 0, the branches a± coalesce and there exists an inte rv al of r over 
which the eigenvalues are complex. S also plays a role in determining parameters in ( ^.12 )-(3.13) other 
than A, e.g. the values (rmid, Cmid) giving the intersection point of the asymptotes of the hyperbola or 
the center of the ellipse. But it is the dependence of A on which leads to the most striking results. 



We now discuss each of (3.12) and (3.13) in turn. Expression ( |3.12 ) can be inverted as follows: 

<J±{l + L) + L 
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■ a±(l + 5) + (5 + L) ^^-^^^ 

(See e.g., [ l|, p[ for similar formulas.) Equation ( [3.15[ ) gives r as a single- valued function of a, despite the 
fact that (3.11) would normally yield an equation quadratic in both a and r. Geometrically, this can be 
understood as follows (see Appendix The roots of the two factors on the left hand side of (3.12) are 
the asymptotes of the hyperbola whose slopes in the (r, a) plane are the eigenvalues Ai± of the matrix 
Ml in (3.7) with determinant Deti. Here, Ai+ = 1 -I- 5, Ai_ = 0, and Deti = 0. Physically, Ai_ and 
Deti vanish because the velocity w advects both the conductive temperature and concentration profiles 
in the same way. Hence a — — {S + L)/{1 + S) is a horizontal asymptote. A line in the (r, a) plane which 
is parallel but not equal to an asymptote intersects the hyperbola in exactly one point (see Appendix 
The consequence of this is that all real values (except — (5* -|- L)/(l -I- S)) of a are achieved exactly once. 
(This distinctive feature of the eigenvalues arising from this problem will also have repercussions on the 
nonlinear problem, discussed in section^.) In particular (except when S = —L), a = for a unique value 
of r, 

r.. = ^ (3.16) 

which is the location of the unique steady bifurcation from the basic state, well-known to be a pitchfork 
bifurcation. For S > (<) — L, the steady bifurcation occurs at positive (negative) r, with rpp +{—)oo 
as 5 i (T) - L. 

If > —L^ or if < — L, then it is the upper branch a+ which crosses zero at the bifurcation, i.e. rpp 
satisfies / -I- ^ = (see cases S = —0.01, —0.1,-1 of figure ||). If S is in the range —L < S < —L^, then 
cr_ crosses zero, i.e. rpp satisfies / — = (see case S — —0.001 of figure At the endpoint S = —L 
of this range, the steady bifurcation goes to infinity, while at the other endpoint — i^, it coalesces with 
a Hopf bifurcation (see below). 



We now consider the ellipse ( |3.13D describing uj{r). If S" < 0, the eigenvalues are complex over the 
range: 

1-L 1-L , , 

r_ < r < r+ = . ^ (3.17) 



(1 + ^/^)2 ^ (1 - V^)2 

The endpoints r± are solutions to g — 0. In this range, a± is equal to half the trace of the matrix M: 

1 



a+=a- = -[r-l + Sir- L)] 



(3.18) 



A Ho pf bifu rcati on occurs at r = Th if cr± = within the range ( 3.17 ) of complex eigenvalues. Equations 
(3.18) and (3.13) show that Th and t^n = '^(''h) satisfy: 



1 + L 
1 + S 



" 1 + S 



(3.19) 
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Thus, a Hopf bifurcation occurs if and only if Th falls within the range (3.17), i.e. 

-1<S < -L^ (3.20) 

If S ~ —1, then Th, r^, and o^h all become infinite. This is the exceptional case in which r{a) forms a 
leftward-opening parabola and r{uj) a rightward-opening parabola (see figure |2|). If: 

S^S^ = -L^ (3.21) 

tl^en: Th = r+ = rpp = = (3.22) 

and io-ti vanishes. This is the well-known codimension-two point, e.g. ^, ||, ||, ^ see figure 

I 

A detailed representation of the bifurcations undergone by the thermosolutal system, made possible 
by logarithmic scaling, is given in figures ^ (for negative S) and ^ (for positive S). 

3.3 Linear thermal and solutal regimes 

From figures || and ||, it can be seen that the thermosolutal thresholds and eigenvalues are related to 
the pure thermal and solutal thresholds and eigenvalues. Such a resemblance is also clearly visible in 
a numerical study of Marangoni (surface-tension-driven) convection with Soret effect ||59|| . This is the 
relationship we wish to explore in this section. 

The simplest classification is by proximity: a real eigenvalue is primarily thermal if it is closer to the 
pure thermal eigenvalue than to the pure solutal eigenvalue CTc, i.e. if: 

|cr - (TtI < |cr - (Tel (3.23) 

Another possible classification is based on the eigenvector, specifically on t he m agnitude of SC/T - the 
ratio of the solutal to the thermal contribut ion in the buoyancy force in ( |3.1c| ) . These two criteria are 
equivalent. Indeed, the eigenvalue equation ( |3.6| ) states that the eigenvectors satisfy: 

SC 

fj — CTt — (3.24a) 
<T-(7c = (3.24b) 



Thus ( |3.23 ) becomes: 

SC 



y <VI^I (3.25) 

Thus, an eigenvalue is thermal (solutal) if the corresponding eigenvector satisfies: \SC/T\ < {>)^y\S\. 

For positive S, rjnt separates thermal from solutal portions of the eigenvalue curve: For < 5 < 1, 
(T+ is thermal and is solutal for the range r > and vice versa r < nnt- For S* > 1, the opposite 
holds. For negative S, rint is the midpoint of the interval of complex eigenvalues. Complex eigenvalues, 
whose real part is equidistant between the pure thermal and solutal eigenvalues, cannot be classified in 
this way. Instead, it is r± which serve as boundaries: for negative S, we classify the real eigenvalues (7+ 
as thermal and ct_ as solutal for the range r > r+ and vice versa for r < r_ . In figures I and |, we can 
discern various regimes in which the thermosolutal eigenvalues a± adhere closely to the pure solutal and 
thermal eigenvalues CTt and ac- For \S\ ^ 1, for example, cr+ w fic for r sufficiently small. For \S\ ^ L, 
a± resemble Ut and ctc over most of the range shown surrounding ri^t- Detailed estimates, omitted here, 
justify these visual impressions. 

Similarity between thermosolutal and pure thermal or solutal eigenvalues implies similarity between 
the thresholds of these problems. The threshold of the steady bifurcation is rpp — L/{L + S). If \S\ <IC L, 
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then Tpp w 1 = Tt, i.e. the thermosolutal threshold approaches the pure thermal threshold. Conversely, if 
\S\ 3> L, then rpp L/ S — rc, so that the thermosolutal threshold approaches the pure solutal threshold. 



More precisely, we can calculate the difference between the thresholds: 

L 



L + S 



- 1 



L-{L + S) 




S 


L + S 




L + S 



1"PF 1"c 1 — 


L L 




LS-L{L + S) 




L^ 




L + S S 




S{L + S) 




S{L + S) 



Taking into account the different sign possibilities for S and L + 5, we find that, for e <C 1, 

|rpp r-p I < e if IS*! < Le 
|rpp-rc|<e if \S\ > ^ 



(3.26a) 
(3.26b) 

(3.27a) 
(3.27b) 



Indeed, on the left diagrams of figures || and |[ we see that the thermosolutal threshold rpp (heavy solid 
curve) is well approximated by the pure thermal threshold (straight thin dashed curve) for \S\ < 0.001 
and by the pure solutal threshold r g (curved thin dashed curve) for \S\ > 0.03. This is precisely the 
estimate obtained from ( 3.27a - 3.27b ) with e = 0.1 and L — 0.01. The domain \S\ < O.IL is analogous to 
what was termed the thermal-dominated regime by Bergeon et al. while the \S\ > 3L is analogous 
to what was termed the solutal-dominated regime. We note further that the eigenvector at the steady 
bifurcation satisfies 

^ - f (3.28) 

Thus, at the bifurcation, either the solutal contribution SC or the thermal contribution T dominates the 
buyoancy force, according on whether one is in the solutal-dominated or thermal-dominated regime. 
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Figure 2: Behavior of thermosolutal eigenvalues for negative S. 

Left: Thresholds for negative values of S plotted on a logarithmic scale. Solid curves show the thresholds 
= L/{L + S) of steady bifurcations. Heavy long-dashed curve indicates the thresholds — (1 + 
L) I (1 + 5') of Hopf bifurcations. This curve appears from r = ooatS' = — 1 and disappears by meeting the 
steady bifurcation curve in a codimension-two point (Bogdanov bifurcation) at = —L^, = 1/(1 — L), 
indicated by a heavy dot. Between the two dotted curves r = r±, the eigenvalues are complex; the right 
boundary r = r_|_ goes to infinity at S = —1. Light dashed curves show the thresholds = 1 and 
Tc — L/S of the pure thermal and solutal problems, which are within 0.1 of rpp for S > — O.IL = 0.001 
and S < — L/VOJ" = —0.03, respectively. 

Right: Real part a± of the eigenvalues of the thermosolutal problem as a function of r for representative 
negative values of S. Straight segments show the real part of complex conjugate pairs of eigenvalues. 
Dots represent bifurcations. Dashed lines are the eigenvalues CTt and ctc of the pure thermal and solutal 
problems, with slopes 1 and S, respectively and intersection point rjnt- For S — —3, representing S < — 1, 
a steady bifurcation with (7+ = occurs at rpp = —0.00334 « —0.00333 = re- Complex eigenvalues all 
have negative real part. For the limiting case S = —I, the range of complex eigenvalues extends from 
r_ = (1 — -L)/4 to r_|_ — CO. All have cr = f — —(1 + L), with uj = \J—g ^ oo as r ^ oo. For 
S = —0.1, representing — 1 < 5 < — L, there is both a Hopf bifurcation at Th ~ 1.12 and a steady 
bifurcation at rpp — —0.11. For the limiting case S = — L = —0.01 (not shown), the pitchfork bifurcation 
has disappeared to +oo. For S — —0.001, representing — L < S < —L^, the pitchfork bifurcation has 
reappeared from —oo, with ct- = at rpp = 1.11. The limiting case S = —L^ = —0.0001 is the 
codimension-two point. Note change of scale between upper two and lower three diagrams. 



4 Nonlinear analysis 

4.1 Derivation of nonlinear model 

We return to the full nonlinear thermosolutal problem ( ^.4| ) and summarize the derivation of a minimal 
set of amplitude equations. This model was first introduced by Veronis [Q for the thermosolutal problem 
and later adapted by Flatten and Chavepeyer H for the Soret problem. 



The simple spatial dependence (3.2) is not preserved by the nonlinear terms in the governing equations 
( ^ . Substituting (HJa) and 

(t){x, z, t) = (t){t) sm{kx) sin(7rz) (4-1) 



where w = k(j>, into (2.4a) yields 



ey(V(/.xVT) = d,(t>d^,T - d^4>d,T 

= —Trkcfyr {sit? [kx) + cos^{kx)) sin(7rz) cos(7rz) 

= -|wrsin(27rz) (4.2) 
and similarly for C . On the other hand, the dependence ( [1.1| ) is preserved by the nonlinear terms of 

q-^ey ■ {\74> X VV^(^) = -Gy • (V0 X V0) (4.3) 
Expansion ( |3.2| ) is therefore generalized to include terms of the type ([1.2|): 

f = T{t)cos{kx)sm{TTz) +T2{t)sm{2TTz) (4.4a) 

C = C(i)cos(fca;)sin(7rz) -|-C2(i)sin(27rz) (4.4b) 
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Figure 3: Behavior of thermosolutal eigenvalues for positive S. 

Left: thresholds for positive values of S plotted on a logarithmic scale. Solid curve shows the threshold 
Tpp = L/{L + S) of steady bifurcations. Dashed curves show the thresholds Tt = 1 and rc = L/S 
of the pure thermal and solutal problems, which are within 0.1 of rpp for S < O.IL = 0.001 and 
S > L/VoJ = 0.03, respectively. 

Right: Real part a± of the eigenvalues of the thermosolutal problem as a function of r for representative 
positive values of S. Dashed lines are the eigenvalues CTt and ctc of the pure thermal and solutal 
problems, with slopes 1 and S, respectively and intersection point rint- Dots represent bifurcations. 
S = 0.0003 illustrates an extreme case of avoided crossing: cr+ is very close to ctc for r < nnt and to CTt 
for r > Tint . This is in the thermal-dominated regime: rpp = 0.97 w 1 = Tt- As S* increases, a± deviate 
more from the pure eigenvalues. For S = 0.03, rpp = 0.25, as compared to Tt = 1 and rc = 0.33. For 
S = 0.3, rpp = 0.032 w 0.033 = rc- For 5 = 3, the slope of CTc exceeds that of CTt. The thermosolutal 
eigenvalues deviate substantially from the pure eigenvalues, but (t_|_ i=a (Tc near the bifurcation point. 
Note change of scale between upper three and lower diagrams. 
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The nonlinear term of (2.4a) on and T2 yields 



3y • (V0 X VT2) = d,^d.,T2 - d^<t>d,T2 
= — 2tiwT2 cos{kx) sin(7rz) cos(27rz) 
= TTwT2 cos(fca;)(sin(7rz) — sin(37rz)) 



(4.5) 



The expansion ( 4.4a - 4.4b ) is truncated by neglecting the term cos(fca;) sin(37rz). This procedure is quan- 
titatively correct for small amplitudes and, as is often the case, qualitatively accurate even for moderate 
amplitudes. The validity and limitations of this truncation are discussed in ||, 



Substituting (H.4a[|4.4bD , ([4.lD , and (|3.2|c) into (|2.l| ) yields: 



d 
di 



( T \ 

C 
w 
T2 

V C2 J 



( -1 



Pr 







PSr 





1 
1 

-P 









-4^2/(72 










-4^2^/(72 J 



( T \ 

c 

w 
T2 

V C2 ) 



T2 
C2 


-r/2 

V -^/2 J 



(4.6) 



In the linear part of ( |4.6| ), (T, C, w) is decoupled from {T2, C2), and the eigenvalues of the latter system 
are always negative. Hence, the linear stability analysis of (4.6) is that already carried out in section ||, 
assuming again that P 00. 



The remainder of section ^ is devoted to analyzing the exact steady states of (f.E). Surprisingly, we 
will find that calculating these steady states reduces to calculating the eigenvalues of a 2 x 2 matrix 
closely related to the linear stability matrix (3.6). 



4.2 Steady states of the nonlinear model 



Although ( j4.6D is a five-dimensional nonlinear system, its form allows its steady solutions to be calculated 
analytically Steady solutions of (4^) satisfy: 



/ ^ 



















V ) 





( 


-1 





1 





\ 


( ^ \ 







-L 


1 








C 




Pr 


PSr 


-P 








w 













-47r2/g2 





T2 


V 














-47r2L/g2 j 


\ C2 I 



TTW 



T2 
C2 



-r/2 

V -^/2 J 



(4Ja) 
(17b) 
(4^c) 
(17i) 
(&) 



Note that (4/7c;), which was derived as (3^) from the time-evolution equatio n un der the assumption of 
large P, is here merely a consequence of the search for steady states and of (O). Using (4/7c;), (4.7d), 
and ([4.7| e) to ehminate w, T2, and C2 yields: 



r — 1 Sr 
r Sr ^ L 



T 
C 



9 2 



1 
i 



T 
C 



(4.8) 



Equation (18) is of the form of an eigenvalue problem: 



r- 1 
Lr 



Sr 

L{Sr~L) 



T 
C 



with 



? 2 



2 2 
q r 



= E 



(T -f SCf 



T 
C 



(4.9) 



(4.10) 



playing the role of an eigenvalue. After the eigenvalues E of ( |4.9| ) are found, the relative amplitude of 
components T and C is given by the eigenvectors. Note that the original five-dimensional steady-state 
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problem (4.7) is also of this type. Then, ~t:w plays the role of a generalized eigenvalue and the generalized 
eigenvectors give the relative proportions of the five components. 

A related version of this reduction has been shown to be valid for the full PDEs governing binary 
fluid convection with Soret effect with realistic boundary conditions, and has been used [ p7| to derive 
a sophisticated and realistic 11-mode model. HoUinger et al. |5^ demonstrate that the velocity 
field can be very well approximated by a single spatial mode and can be adiabatically eliminated. The 
nonlinearities are then only those which advect the temperature and concentration, and contain only the 
amplitude w of the vertical velocity. The temperature and concentration amplitudes are then solutions 
to a linear system depending on w. HoUinger et al.'s work provides evidence that our reformulation 



of (4.7) as an eigenvalue problem, while strictly valid only for the minimal five-mode Veronis model of 



thermosolutal convection, is a manifestation of a quite general property of double-diffusive convection. 
E can plausibly be called the energy of the convective state, since the kinetic energy density is: 



1 



k f'^'^^^ 

— I dx I dz— [ ( — ^ sin(/ca;) cos(7rz) ) -f (w cos(fca;) sin(7rz))^ 



J 
8fc2 /c2 



2 V V fc 



(4.11) 



The most common experimentally measured quantity is the convective heat transport A/" — 1, which is 
related to E via: 



AA-1 



27r/fc pi 

dx dz w T 
Jo 

2 



TT TT W 
W 1 ^ — 

2k 2kE+l kq^ E+1 



2iv/k j.1 

dx dz w T cos^(fca;) sin^(7rz) 
Jo 

47r E 



(4.12) 



where we have used w — T ^ ET from (4.9) 
The matrix: 



r — 1 Sr 
Lr L{Sr - L) 



M 









-i2 



Afo + rMi 



1 S 
L LS 



(4.13) 



is very similar to the matrix (3.7). Mq again describes diffusion, and Mi now describes the combined 
effects of advection and nonlinear saturation. The interpretation of the behavior of the eigenvalues of 
(4.13) is, of course, different from that of (3.7). In the linear system, negative or complex eigenvalues 
characterize infinitesmal perturbations which decay and/or oscillate. For the nonlinear problem, negative 
or complex values of E are forbidden by virtue of definition ( 4.10 ), and imply non-existence of steady 
solutions for certain ranges of r, S, L. Since only real positive values of E are meaningful, we will not 
introduce separate notation for real and imaginary parts of E. For E real and positive, we define 
A — ±Ve as the convection amplitude. 



4.3 Bifurcations 

Most of the results of section || concerning the eigenvalues a± are easily modified to apply to the energy 
E of the nonlinear steady states, merely by substituting 

L (4.14a) 
S^LS (4.14b) 
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r r 



Figure 4: E and A = iV^E (for E real and positive) for S = —3 and S = 1. For S = —3, E undergoes 
complex coalescence. The point at which E becomes complex corresponds to a pair of saddle-node 
bifurcations in A. For S — 1, E undergoes avoided crossing. This leads to a distinction in A between 
low-amplitude Soret and high-amplitude Rayleigh regimes. Pitchfork bifurcations are at r = 0.0099 for 
5 = 1 and at r = -.003344 for S = -3. 



When we use the same notation for the nonlinear and linear problems, we add tildes to designate nonlinear 
quantities. 

Just as we did for the linear problem, we can define nonlinear pure thermal and solutal solutions by 
setting the coupling terms in (O) to zero: 

ET = r- l= r-rT (4.15a) 
Ec = L{Sr - L) = LSir - rc) (4.15b) 

Again, the coupling term 

(3j = LSr"^ (4.15c) 

is proportional to, and has the same sign as, S. Hence the eigenvalues E undergo avoided crossing near 
the intersection point of the two pure solutions 

if S is positive, and complex coalescence if S is negative. These possibilities are illustrated in figure |[ 
along with the resulting consequences for A = ±\/E. 



15 



The solutions to the eigenvalue problem (4.9) are: 



- . ,,,, 

= \[{l + LS)r-{l + L^)] 

±i V(l + LSfr'^ - 2(1 - i2)(l - LS)r + (1 - L2)2 



/(r) ± 



(4.17a) 

(4.17b) 
(4.17c) 



This expression for E is quite similar to ( B.llD for the eigenvalues. Detailed results for representative 
values of S are shown in the two large figures |5| and o. 



Just as A = 0, (T, C) = (0, 0) is a solution of (|3J), it is also true that E = 0, (T, C) = (0, 0) is always 
a solution to (4.9). The existence of the solution i? = for all r reflects the fact that the conductive 
profile in a motionless fluid rema ins a solution for all Rayleigh numbers, although not necessarily a stable 
one. Equations ( 3.11 ) and ( 4.17 ) give the nontrivial solutions for A and i?, which are zero only at isolated 
values of r. 



The values ( 4.17 ) for E real form hyperbolas: 



E 



LS + L^ 
1 + LS 



-{l + LS)[r- 



{l-L')il-LS) 
{1 + LSy^ 



E 



LS 



1 + LS 



(4.18) 



in the (r, E) plane. (In the exceptional case S = — 1/L, the solutions form a parabola.) Just as we found 
that r was a single- valued function ( |3.15 ) of a, here r is also a single- valued function of E. (This is again 
because the zero determinant of matrix Mi in ( 4.13| ) leads to a horizontal asymptote for the hyperbola.) 
We use the substitutions (4.14) to transform r{a) into r{E): 



<y± + <7±{i + L) + L 

a±{l + S) + {S + L) 



E^ 



E{l + L^) + L^ 



E{1 + LS) + L{S + L) 



(4.19) 



a result which can verified by inverting ( [1.17] ). Similar formulas are found in, e.g. ^, |T^. Thus each 
value of E except —L{S + L)/{1 + LS) is achieved exactly once, a feature mentioned in section |^ for the 
growth rates a, but more important in the context of amplitudes of steady states; see figures I and |. It 
is remarkable that this property continues to hold even for the far more complicated Soret problem with 
rigid boundaries and adequate spatial resolution |5^. In |57|, it is also found that r is a single- valued 
function of the energy, as measured by the square of the vertical velocity amplitude. This function, like 
( 4.19 ), is the ratio of a quadratic to a linear function of E, although with coefficients more complicated 
than those of (|4.19| ). 



Steady bifurcations from the conductive profile occu r at v alues of r at whic h E = 0. Equation (3.16) 
gives the values rpp at which cr = 0. The substitutions ( 4.14 ) leave expression ( 3.16| ) unchanged: 



L 



L2 



L 



S + L 



rpF 



LS + L^ S + L 



(4.20) 



This result can be verified by setting E = in ( 4.17 ) or ( 4.19 ). The fact that the "linear" (cr = 0) and 
"nonlinear" {E = 0) values of Tpf are identical reflects the fundamental bifurcation-theoretic fact that 
a change of sign of a real eigenvalue signals the transverse intersection of two or more steady branches, 
i.e. a steady bifurcation. The fact that E — (rather than A) changes sign at rpp identifies these as 



pitchfork bifurcations, since either sign of A is permitted. Equations (4.12) and (LIE) can be used to 
calculate the slope of the convective heat transport at the bifurcation: 



d(AA- 1) 



dr 



47r 1 

V (£; + 1)2 



dr 
dE 



E=0 



Att L{S + L)^ 
W S + L3 



(4.21) 
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very similar to the analogous formula derived for the Soret problem in |pl[ 



If 5 < 0, complex solutions to (4.17) are possible. The interval of r- values over which this is so, 
delimited by g — 0, is easily obtained using substitutions (4.14) in (3.17): 



1-L2 



(1 + V^)2 



(l-^/=I5)2 



(4.22) 



For the linear problem, the endpoints of the complex interval ( 3.17 ) for A divided oscillatory evolution 
from monotonically growing (if a > ) or decaying (if cr < 0) evolution. For the nonlinear problem, 
the endpoints of the complex interval ( 4.22| ) mark the appearance or disappearance of real solutions A; 
see figure I for < and figure | for 5 = -10. These are saddle-node bifurcations: the simultaneous 
creation of four branches of steady solutions A, corresponding to two different values oi E = A^. Since 
E must be non-negative as well as real, the occurrence of saddle-node bifurcations at f_ or f+ requires 
that n{E) = / > 0. Substituting (f4.22| ) into 



/>±) 



1 



we calculate that saddle-node bifurcations occur at Vs 

1 



= r+ if 



(4.23) 



(4.24) 



Figure |^ characterizes the saddle-node bifurcation over the domain —1/L<S< —L^ of its existence. 
Over most of the range, the value E^j^ = E{rsN) varies with S like l/\/—SL and varies like l + 2\/—SL. 
When the saddle-node bifurcation appears at 5' = — 1/L, we have rg^ = E'gn = oo. When the saddle- node 
bifurcation disappears at the degenerate pitchfork at 5 = —L^, we have — 1/(1 — -^^)j -E'en = 0. The 

■2, 



curvature d r/dA (rsfi) at the saddle-node bifurcation is also shown, d r/dA (rs^) ~ 4 over much of 
the range. In contrast, the curvature at the pitchfork bifurcation d'^v/dA'^i^Vpp^ ^ (^di^Af — l)/dr) "'^(rpp) 
(see equations ( 4.19| ) and ( 4.21 )) varies over many orders of magnitude (even away from its divergence 
at 5 = —L). This great difference in curvatures is one of the factors giving the thermosolutal bifurcation 
diagrams, such as the S = —0.003 case in figure |^, their characteristic appearance. 



We note that ( 4.24 ) is the transformation via substitutions ( 4.14 ) of the range of existence ( 3.20 ) 
of Hopf bifurcations of the linear problem. Yet the criterion for a Hopf bifurcation, where A has zero 
real part (/ = 0) and finite imaginary part {g < 0) has no significance for a nonlinear steady state E. 
We explain the correspondence between the range of existence of Hopf and saddle-node bifurcations as 
follows. 

Although a Hopf bifurcation has no analogue for nonlinear steady states, the codimension-two point of 
the linear analysis, at which the pitchfork and Hopf bifurcations coalesce, does have a relevant nonlinear 
analogue. This is because at the codimension-two point, the pitchfork and Hopf bifurcations necessarily 
coalesce with a third point: that at which the eigenvalues become complex. (See S — case of figure 
^.) That is, the conditions defining a steady bifurcation, / ± y/g = 0, g > and a Hopf bifurcation, 
/ = 0; 5 5: together imply f = g = 0. For the nonlinear analysis, the codimension-two point marks the 
coalescence of the pitchfork bifurcation, defined by / ± ^/g — 0, g > 0, and the saddle-node bifurcation, 
defined by / > 0, g = 0. These conditions together again imply f = g = 0. 



We use (4.14) to transform the codimension-two point of the linear problem to the codimension-two 
point of the nonlinear problem. 



1 



LS^ = -L^ i.e. 
1 



- — - L 



1 -L 



r* = 



1 -L2 



(4.25) 
(4.26) 



Whereas the linear codimension-two point 5* = —L^ marks the limit of existence of the Hopf bifurcation, 
the nonlinear codimension-two point 5"* ~ —L^ marks the limit of existence of the saddle-node bifurca- 
tions. For S > —L^, the pitchfork bifurcation is forwards, with a pair of nontrivial solutions branching 
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right towards r > rpp. For S < —L^, the pitchfork bifurcation is backwards, with nontrivial solutions 
branching left towards r < rpp. 



The correspondence between the lower endpoint of ranges (3.2C) and (4.24) is also easily explained. 
As S I —1, the Hopf bifurcation point Th diverges to +00, as seen in figure |^. This event coincides with 
the divergence of the right endpoint ?■+ of the interval of complex eigenvalues A, as the set {r,±uj(r)) 
evolves from an ellipse to a righward-opening parabola. Similarly, the right endpoint f+ = of the 
interval of complex (forbidden) values of E diverges to +00 as LS i — 1 or | — 1/L, as seen in figure ^ 

The saddle-node bifurcation for the nonlinear thermosolutal and the Soret problems is well-known, 
as are the Hopf bifurcation and the codimension-two point of the linear problem. Indeed, these are the 
basic features that originally inspired the great interest evoked by binary fluid convection. However, 
the relationship between these phenomena has not been previously formulated. One of the advantages 
of the idealized free-slip thermosolutal model is that the codimension-two points have the simple forms 
5"* = — and 5* = —L^. However this scaling is quite general. For the Soret problem with free-slip 



permeable boundaries and finite Prandtl number P ||7|, |2l], |2^ , 

^* = P + L{f+L){P+l) ^* = 'l + L^L^ + L^ ^^-^^^ 

Schopf and Zimmermann |50| observe the scaling S^, = —L^ and S^, = —L^ in amplitude equations 
calculated for the Soret problem with rigid impermeable boundaries. Using their 11-mode minimal 



model, HoUinger et al. |57 calculate 



,2 0.35P-I-0.18 1.96 



^" P + 0.28L ^' 1 + 1.04 L + 0.37 L2+ 0.97 ^^■'^^^ 

Exact expressions for the lower bounds for the existence of the Hopf (5* > — 1 in our case) and saddle- node 
{S > — 1/i in our case) bifurcations are not as readily available, but also continue to hold approximately. 
This provides empirical evidence that the relationship between the complex eigenvalues of the linear 
growth-rate problem and the saddle-node bifurcation of the nonlinear steady-state problem continues to 
hold, at least approximately, even for a more complicated and realistic case, and thus that this relationship 
is a fundamental feature of binary fluid convection. 
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Figure 5: Bifurcation diagrams for negative S. 

Left: thresholds for negative values of S plotted on a logarithmic scale. Solid curves show the thresholds 
= L / {L + S) of steady bifurcations. Dotted curve shows the thresholds rsN for saddle- node bifurca- 
tions. This curve appears from r = +00 at S = —1/ L and disappears by meeting the steady bifurcation 
curve in a codimension-two point (degenerate pitchfork) at 5* — ~L^, — 1/(1 — L^), indicated by a 
heavy dot. Dashed curves show the pitchfork bifurcation thresholds r^. — 1 and = L/{L + S) of the 
pure thermal and solutal problems. 

Right: Energy E and amplitude A of steady nonlinear solutions to thermosolutal problem as a function 
of r for representative negative values of S. Dashed curves on lower two sets of diagrams show energy 
and amplitude of pure thermal and solutal solutions. For S — —300, representing S < —1/L, there is one 
branch of real solutions A, bifurcating towards negative r and resembling the pure solutal branch. At 
S = — 1/L (not shown) , the set (r, E) is a parabola and a pair of saddle-node bifurcations descends from 
r — 00. For S = —10, representing —1/L < S < —L, the resulting pairs of disconnected branches can be 
seen inside and rather close to the pure thermal branch. The other solutions, branching towards negative 
r, are still present. They, and the pure solutal solutions which they resemble, are of greatly decreased 
amplitude. At S* = —L (not shown), the pitchfork bifurcation disappears at r = —00 to reappear at 
r — +00. (The pure solutal branch, increasingly small and distant, continues to exist until S = 0.) For 
S = —0.003, representing —L < S < —L^, the pitchfork bifurcation has descended to r = 1.429 and 
connects the two pairs of branches arising from the saddle-node bifurcations. Here, and in the next case, 
the thermal branch resembles the upper branch of solutions too closely to be distinguished from it in the 
figure, whereas the solutal branch cannot be distinguished from the r-axis. At 5' = —L^ — — 10~^, the 
saddle-nodes and pitchfork coalesce in a codimension-two point separating subcritical from supercritical 
pitchfork bifurcations. Note change of scale between upper two and lower two diagrams. 



4.4 Soret and Rayleigh regimes 



In section 3.3, we saw that the eigenvalues could be characterized as primarily thermal or primarily solutal 
according to their distance from the pure thermal or solutal eigenvalues. Here we will discuss various 
ways of classifying the nonlinear solution branches in this way. This classification is more significant 
since we are interested in entire branches of nonlinear steady states, whereas eigenvalues are of interest 
primarily at the thresholds. 

We can apply the classification by proximity analogous to that we used for the eigenvalues: a nonlinear 
steady state is primarily thermal if it is closer to the pure thermal than to the solutal branch, i.e. if 



\E-Et\ < \E^Ec\ 



(4.29) 



and primarily solutal otherwise. Just as we did for the eigenvalues and eigenvectors, we can show that 
this criterion is equivalent to one based on the magnitude of SC/T, which is the ratio of the solutal to 
the thermal contribution in the definition of the convective amplitude A cx T -I- SC as well as in the 
buoyancy force. The eigenvalue equation (K^) states that the nonlinear steady states satisfy: 



E- El 



SC 
T ^ 
LT 

E-Ec = —r 



Thus ( [4.29D becomes: 



SC 



T 



< 



^/\LS\ 



(4.30a) 
(4.30b) 

(4.31) 



and so a nonlinear steady state is thermal (solutal) if SC/T < {>)^/LS. 

We first consider negative S. Four ranges of S can be distinguished, as can be seen on figure ^. For 
—L^ < S < 0, a pair of branches bifurcates towards positive r; these nonlinear steady states are all 
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Figure 6: Bifurcation diagrams for positive S. 

Left: Thresholds for positive values of S plotted on a logarithmic scale. Solid curve show the thresholds 
= L/{L + S) of steady bifurcations. Long-dashed curve shows value fint at which the pure thermal 
and solutal branches intersect. Dotted curve shows value fmid at which the asymptotes of the hyperbola 
E intersect, and at which the slope of E and the curvature of A change most rapidly. 
Right: Energy E and amplitude A of steady nonlinear solutions to thermosolutal problem as a function 
of r for representative positive values of S. Dashed curves show energy and amplitude of pure thermal 
and solutal solutions. At S' = 0.001, representing < S < L , the thermosolutal branch hugs the 
thermal branch. At 5* = 1 and S = 10, representing L < S < 1/L, the pure solutal branch has descended 
below the pure thermal branch, "pulling" the thermosolutal branch with it. At S = 100, represent- 
ing S > 1/L, the solutal branch lies above the thermal branch. The thermosolutal branch, initially 
tangent to the solutal branch, lies above both pure branches. All bifurcation diagrams use the same scale. 
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-L-1 -1 -U -L2 -L3 
S 



Figure 7: Saddle-node bifurcation for —1/L < S < ~L^. The solid curve shows £'sn, the energy at 
the saddle-node bifurcation at Tsn- -Esn ~ \/\LS\ except very near — where it diverges, and near 
—L^, where it goes to zero. The dotted curve shows (Pr/dA^lrsn), the curvature of the saddle-node 
bifurcation, which remains near 4 except near the endpoints —1/L and —L^. For contrast, the dashed 
curve shows (Pr/dA'^{rpF), the curvature of the pitchfork bifurcation, which varies greatly everywhere 
over this interval, especially near — L, where rpp and cPr/dA'^{rpp) approach ±cxd. 



thermal. For —L < S < —L^, the pitchfork bifurcation is backwards and both and E- are real and 
positive for < r < rpp. Thermal and solutal branches are separated by f+ — Tsn, where SC/T = \J LS. 
The lower branch E^ is solutal (but there exist no corresponding pure solutal steady states, since Eq is 
negative) and the upper branch E+ is thermal. For — 1/i < S < —L, the pitchfork bifurcation occurs 
at negative rpp and branches towards lower r; the resulting branches are solutal. (These do correspond 
to pure solutal steady states.) The branches that exist for r > Tsn are isolated. The high-amplitude 
branches £'+ are primarily thermal and the low-amplitude branches i?_ solutal. For S < —1/L, the 
isolated branches no longer exist, leaving only the solutal branches, whose amplitude increases with \S\. 

We now focus on positive S. Although less studied than negative S, this case has nonetheless re- 
ceived substantial attention for the Soret problem. Qualitatively, in experiments or three-dimensional 
calculations, a striking pattern of squares is produced (2|,|o[||||,|o[|5[||,|4[|||o[|ll, possibly 
alternating with rolls of different orientation Quantitatively, a fairly abrupt transition is observed be- 
tween a low-amplitude and a high-amplitude convcctive regime, as can be seen in figures ^. This abrupt 
transition was first derived for the five- variable model by Flatten and Chavepeyer first observed 
experimentally by Le Gal et al. ||2^, and the two regimes identified and named the Soret and Rayleigh 
regimes by Moses and Steinberg |3C|, The Soret-to- Rayleigh transition has also been reproduced 

experimentally in [^9[ |4| and numerically in |4j, ^ |5^. We will interpret this transition as a 
manifestation of the avoided crossing phenomenon which occurs at the separation between the solutal 
and thermal regimes. 

The physically significant real and positive values of E are those for r > rpp- Thus, a steady state 
branch has physically significant solutal and thermal portions if fjnt > ''pf, which occurs if < S < 1/L. 
In this case, the solutal regime comprises: 

^ = rpp < r < nn, ^ 1-^ (4.32) 
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The transition between solutal and thermal steady states is manifested by the related increase in slope 
(for E) or curvature (for A). For S > 1/L, the entire solution branch is solutal while ioi < S < the 
entire solution branch is thermal (see figure |^) . 

What is actually observed is more complicated than the analysis given above. Recall that the asymp- 



totes of the hyperbolas (4.18) describing E differ from the lines Exir) and Ec{r). The increase in slope 
of E occurs, not at the intersection point fi„t between the pure solutal and thermal branches, but at the 
intersection point fmid between the two asymptotes; see Appendix For < 5* < 1 and for S > 10^, 
fmid and fint are so close as to be indistinguishable on figure ^. The same proviso holds for f-aiid as for 
fint^ the transition is observed only if it occurs at a real and positive value of E, i.e. only if fi^id > fpF- 
In this case, the Soret regime comprises: 

- rpp < r < r^id = ^ ,^ \ (4.33) 



L + S " (i + Lsy 



For S' = 0, we have fmid = 1 — < 1 = rpp, so the entire solution branch is in the Rayleigh regime. 
For S large, it can also be shown that fmid < rpp; the entire branch is then in the Soret regime. But there 
is an intermediate range of S over which the transition can be observed, as can be seen from setting rpp 
equal to fmid: 



L + S 




{1-L^){1-LS) 

{1 + LSf 
LS'^ - (1 - 4L^ 



5*1,2 — 



2L 



(4.34) 



Equation (4.34) has two real solutions 5*1, 52 if 1 — 6L^ + is positive, which occurs if L < VS - 2V2 = 
0.41 or if L > ^/3 + 2V2 = 2 .41, as illustrated in figure ||. For L satisfying these conditions, and for 
Si < S < S2, there exists a low-amplitude Soret regime over interval (4.33). For L <C 0.41, we have 
Si ^ and 5*2 ~ l/L and vice versa for L ^ 2.41. (Recall that fint > ^pf for < S < 1/L.) Figure 
^ shows the variation with S of fmid, ''pf, and fint for the case L = 0.01. We see that while the size of 
the Soret range fmid — ^pf is positive for i"^ < 5 < 1/L, it is appreciable only over the smaller interval 

<S< 

For the transition to be observed, the energy must also be sufficiently large for convection to be 
detected. 

— -L{S + L) + y/SL{l-L^) 



l + LS 



(4.35) 



is maximal at 5* = (1 — L)^/(4L) w 1/(4L) for L <C 1 and is negative or small outside the range 
LlesssimS < 1/L (see, e.g., figure ^ for S=0.001). 

Another factor that bl urs t he transition from Soret to Rayleigh regime is the fact that as S increases 
above 0, the hyperbola (4.18) separates from its asymptotes, and the change in slope becomes more 
gradual as seen, for example, in figure ^ for S=10. (This is somewhat counterbalanced by the fact that 
the angle between the two asymptotes increases from tt/4 at S — to 7r/2 at S — 00.) The change in 
slope of E is measured by E" , which is maximal at fmid (see Appendix Normalizing -E"(fmid) by 
£'(^mid) and taking the inverse square root defines a length in r over which the change in slope occurs. 
Dividing this length by fmid ~ ''pf yields a ratio which compares this length to the extent of the Soret 
regime. Thus we define 



1 



-,1/2 



where 



^^"(rmid) = 



.^"(rmid). 
{1 + LSf 



4V5l(l - L2) 



(4.36) 



(4.37) 
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Figure 8: Range (5i,iS'2) of values for which f,„id > r-pp 
L < 0.41 or if i > 2.41. 



, as a function of L. This range exists for 
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Figure 9: Quantities determining existence of Soret-to-Raleigh-regime transition for positive S. Upper 
diagram: Soret regime exists for rpp < r < f^id. ?^mid — ''pf is positive for < S < 1/L and appreciable 
for < 5- < L-3/4. Lower diagram: Abruptness of transition is measured by the smallness of the 
ratio S (dashed curve) defined in ( 4.36| ) which remains less than 0.2 for for ^ S < 1. Energy E{f t^id) 
at transition (solid curve) is maximal at 1/(4L) and non- negligible for L < S <1/L. 
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and TpF, fmid, and -E(rmid) are given in (4.33) and (4.35). Figure |9| shows that S is smallest - i.e., 



the change is most abrupt - near S = L and rises steeply for 5 > 1 and S L^. The interplay of 
counterbalancing criteria demonstrates the multiple roles played by S in this geometric analysis. 

Combining all of these criteria, we finally obtain L < 5 < 1 as the separation parameter range for 
the Soret-to-Rayleigh transition to be observable. The experimental observations of the Soret regime 
have indeed been approximately in this range. Moses and Steinberg [E5|, who have carried out the most 
extensive experimental investigation, observe the Soret regime for i'^°°< S < L*^'^^. Other experimental 
observations are at 5 = L^-^^ [|9|, at S* = L^-^^ and at 5 = Q. Numerical observations 

should be possible over a larger range: since r and E ranges can effectively be magnified as required, 
the thresholds for r^id ~ ^pf and E{r^i({) are not as constraining. Flatten and Chavepeyer pO| observe 
the Soret regime for L^-^^ < S < LP-^ , other subsequent investigators reported a transition at S* = i^-^^ 
!§, at S" = L"-2i 111, at S" = L^-^e [g|, and at 5 = L^-^ ||. For aU values S" < 1, the transition point 
^mid is indistinguishable from the thermal threshold tt — ^ (see figure ||) and thus the Soret-to-Rayleigh 
transition is invariably described as coinciding with the onset of thermal convection in a pure fluid. 

All of the references cited above have investigated the Soret problem. There, the no-flux boundary 
conditions on C lead to /c = as yielding the lowest threshold for linear instability for sufficiently large S*, 
e.g. This zero-wavenumber instability is sometimes invoked as part of the explanation for the weak 

heat transport in the Soret regime p] , However, the wavenumber actually realized in full nonlinear 
simulations for the Soret problem with rigid boundaries in a two-dimensional geometry is close to 
TT, as is the wavenumber for the square patterns observed experimentally, e.g. Considerations 
of pattern selection clearly play no role in the mechanism we have discussed since, in the thermosolutal 
problem, the boundary conditions on T and C are identical and we have fixed k at fccrit = 7r/-\/2- 

Finally, we consider the asymptotic behavior of the thermosolutal solution branches. For \r\ large, we 
have, for the upper branch 



SC 
E 

'e^ 

E 



LS with \LS\^y/\LS\ if \S\^1/L 



{l + LS)r _ . 
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if \S\ 
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l/L 



if \S\ 



> 



l/L 



(4.38a) 
(4.38b) 

(4.38c) 



Relati on (|4.38a ) shows that the upper branch is thermal if \S\ < l/L and solutal if jS"! > l/L. Rela- 
tions (4.38b- r.38c| ) suggest another, more stringent, criterion for classification: a steady state is thermal 
(solutal) if the ratio E/Et (E/Ec) is close to one, which is true for the upper branch if \S\ <C (^)1/L. 



For the lower branch i?_ 



— « -1 with 1 ^ ^/\LS\ if 1^1 ^ l/L 



(4.39a) 



so the lower branch is solutal if \S\ < l/L and thermal if |5| > l/L, by the criterion ( 4.29| ) but it fails to 
meet the more stringent criterion since the ratios E/Et and E/Eq both tend to zero: 



E 

E 
Ec 



-L{L + S)/{1 + LS) 
r-1 

-L{L + S)/ jl + LS) 
LS{r - S/L) 



(4.39b) 
(4.39c) 



The asymptotic behavior ( 4.38a ) and ( 4.39a ) is illustrated in figure [T 
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Figure 10: Solid curves show T/C (above) and T/{SC) (below) for S = ±1, ±0.1, ±0.01; long- (short-) 
dashed curves show T/C or T/{SC) for S = ±(— )10. As r becomes large (and along the upper of the 
two branches for S negative), T/C approaches its asymptotic value of 1/i, independent of 5*. For low 
amplitude states existing for S negative, T/{SC) approaches —1. 



5 Time-dependent model 



The linear and nonlinear equations of sections |^ and ^ : 

1 
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di 
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r- 1 
r 
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1/L 
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(5.1a) 
(5.1b) 



each lack an essential feature of the bifurcation diagram for binary fluid convection: the linear stability 



problem (5.1a) cannot contain saddle-node bifurcations, whereas the nonlinear steady-state problem 
( ^.lb| ) cannot describe Hopf bifurcations. We can combine (5.1a) and (5.1b) to form a single time- 
dependent two-variable system containing all of these features: 



d_ 

di 



T 
C 



1 Sr 
Sr-L 



T 
C 



-\{^)\t,scY- 



1 \ / T 
l/L )\C 



(5.2) 



The conditions for the validit y of this t wo-di mensional time-dependent nonlinear system combine those 
required for the two system s (^.la ) and ( 5.11 ): large Prandtl number and small amplitudes. In addition, 
the spatial representation (|3.2| ) imposes a fix ed p hase on the solutions. Beyond these two statements, 
we m ake no claim for the accuracy of system (|5.2| ) as a representation of the partial differential equation 

(111)- 



By construction, system (5.2) undergoes a pitchfork bifurcation at r = rpp — L/{L + S) and a 
Hopf bifurcation at Th = (1 + L)/{\ + S"), both from the trivial state, and saddle-node bifurcations at 
r = TgN = (1~^^)/(1~ ^f—LSy. In addition to reproducing the linear stability of the conductive state 
and the nonlinear steady states, this system also displays an interesting phenomenon that can occur in 
actual binary fluid convection: the limit cycle disappears via a global bifurcation by colliding with the 
saddles on the lower branch of steady states. In figure we show numerically computed phase portraits 
of (U) for S = -L = -0.01. Initial conditions are (T, C) = ±(0, 0.2), ±(0, 2.0). For this value of S, the 
Hopf and saddle- node bifurcations occur simultaneously at r — 1.02. For r slightly less than r — 1.02, 
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Figure 11: Phase portraits illustrating heteroclinic bifurcation for S = —L = —0.01. For r = 1.00, 
all trajectories spiral into {T,C) — (0,0). For this value of S, the Hopf and saddle- node bifurcations 
occur simultaneously at r = 1.02. For r = 1.05, trajectories originating close to zero spiral out to a limit 
cycle, while trajectories originating sufficiently far from zero terminate on one of the stable steady states 
(solid dots), possibly after being deflected by one of the saddle points (hollow dots). For r = 1.06, the 
limit cycle has been destroyed by colliding with the saddle points in a heteroclinic bifurcation and all 
trajectories terminate on one of the stable steady states. 



the origin is a stable spiral node; all trajectories spiral into (0,0). At r — 1.02, a limit cycle and two 
pairs of steady states (stable nodes and unstable saddles) are created. For r slightly exceeding r = 1.02, 
trajectories approach either the limit cycle or the stable steady state, depending on the initial condition. 
At r = tg with 1.05 < tq < 1.06, the heteroclinic bifurcation destroys the limit cycle. For r > 1.06, all 
trajectories terminate on one of the stable steady states. The influence of the saddles on the trajectories 
is clearly seen. 



In figure [T^ we show the thresholds for the bifurcations undergone by model (5.2). The accompanying 
bifurcation diagrams are schematic, unlike those of figures |, |, |, and |. There are six qualitatively 
different diagrams, for S in ranges S < -1/L, -1/L < S < -1, -1 < S < -L, -L < S < -L'^, 
< S < —L^, and —L^ < S. This illustrates the advantage of the simplified model: each change in 
the qualitative dynamics occurs exactly at a power of L. Each steady branch is labeled with the number 
of eigenvectors to which it is unstable, i.e. the number of eigenvalues with positive real part. In section 
^, we characterized the pitchfork bifurcations as forward or backward facing, according to whether the 
new solutions created branch towards r > rpp or r < rpp, rather than as supercritical or subcritical. The 
reason for this is that the criterion we use for super or subcriticality depends on a combination of linear 
and nonlinear information: a bifurcation is supercritical if the new solutions branch in the direction of 
increasing instability of the parent branch. In this sense, the pitchfork bifurcations are supercritical for 
all cases except — < 5 < —L^, despite being backward facing for the five cases S < —L^. For the three 
cases S < —1/L, —1/L < S < —1, and —1 < S < —L, a real eigenvalue becomes and remains positive 
as r is decreased below rpp. For the fourth case —L < S < —L'^, the positive eigenvalue coalesces with 
another positive eigenvalue to form a complex conjugate pair whose real part then reverses direction and 
becomes negative as r is decreased, resulting in the Hopf bifurcation at Th, as was shown in figure ^ for 
S = —0.001. Although we have shown the Hopf bifurcation as supercritical in figure ^ we do not exclude 
the possibility of a subcritical Hopf bifurcation accompanied by a stabilizing saddle-node bifurcation. 

In a two-variable system, the heteroclinic bifurcation is a natural consequence of Hopf bifurcations 
and pitchfork bifurcations which branch towards each other. Starting from the Hopf bifurcation and 
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Figure 12: Behavior of the hybrid model for negative S. 

Left: Thresholds for negative values of S plotted on a logarithmic scale. Solid curves show the thresholds 
= L/{L + S) of pitchfork bifurcations, negative for S < —L and positive for S > —L. Dash-dotted 
curve indicates the thresholds re = (1 + L)/{1 + S) of Hopf bifurcations, appearing from r = oo at 
S = —1 and disappearing by meeting the pitchfork bifurcation curve in a codimension-two (Bogdanov) 
bifurcation at S" = — (large dot). Dotted curve indicates the thresholds Tsn of saddle- node bifurcations, 
appearing from r = oo at 5* = — l/i and disappearing by meeting the pitchfork bifurcation curve in a 
codimension-two (degenerate pitchfork) bifurcation at S = —L^ (large dot). 

Right: Schematic bifurcation diagrams for the six qualitatively different cases. Numbers above branches 
indicate the number of eigenvectors to which the branch is unstable. Stable branches (0 unstable 
eigenvectors) are additionally shown as heavy curves. Dotted curves indicate limit cycles. For S < —1/L, 
the pitchfork bifurcation is supercritical and branches backwards from negative rpp. For — l/i<5<— 1, 
additional disconnected branches arc created via saddle-node bifurcations. For —1 < S < — L, a Hopf bi- 
furcation creates a limit cycle which terminates via a global heteroclinic bifurcation. For —L < S < —L'^, 
the supercritical pitchfork bifurcation branches backwards from positive rpp and connects the branches 
emanating from the saddle-node bifurcations. For —L^ < S < —L^, the Hopf bifurcation no longer exists 
and the pitchfork bifurcation is now subcritical. For —L^ < S, the saddle-node bifurcations no longer 
exist and the pitchfork bifurcation is again supercritical. 
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approaching the pitchfork bifurcation, the amphtude of the hmit cycle increases while that of the steady 
branch (of saddle points) decreases. Confined to a plane, the limit cycle and steady states then collide 
at some intermediate value of r. 

More generally, this bifurcation sequence has been studied extensively |l^, ^ . The 

approach to the heteroclinic bifurcation has been observed in numerical simulations of the full system 
of governing partial differential equations for thermosolutal convection ||Tl|, and for the Marangoni- 



Soret problem |59 . The existence of the heteroclinic bifurcation has been proven using the techniques 
of normal form reduction and amplitude expansions, and analytic expressions calculated for the limit 
cycles and the bifurcation ^ In addition, simulations of the five- mode Veronis model exhibit 
interesting complex dynamical phenomena such as period-doubling and chaos p^ ; this is not possible 
for a two- variable model such as (5.2). Mathematical analyses have rigorously derived reduced models 



of binary fluid convection, analyzed their domains of validity, and determined when period-doubling and 
chaos may occur WA B . 



Referri ng t o (3.2), the limit cycl e o f figure is a standing wave solution of the partial differential 
equations (2.1). In treatments of (3.2) adapted to large or infinite horizontal domains, in which the 
phases of the various components are allowed to vary, then the Hopf bifurcation gives rise to a branch 
of stable traveling waves which disappears via a drift bifurcation by meeting the branch of stable steady 
states; the standing waves described above continue to exist but are unstable jl^, |2^, HJ, |7j. Traveling 
waves are indeed observed in containers which are large or periodic in the horizontal direction [p5| 
p8| , |36| , [5^ . The situation is in fact far more complicated: the traveling or standing wave branches may 
bifurcate subcritically, e.g. 1 1 1 , [l5 , ^ ^ , the traveling wave branch undergoes a secondary bifurcation to 

13l, and, at least for the Soret problem, the traveling wave branch can 

" 55! . As a further complication, 
" PI. Yet, 



modulated traveling waves |27, |36, 

undergo several saddle-node bifurcations between slow and fast branches 
in larger two-dimensional domains, localized traveling waves and pulses predominate 
certain large-scale aspects of steady-state convection examined in section ^ have counterparts for standing 
and traveling waves. Schopf and Zimmerman have found that the degenerate Hopf bifurcations for 
the standing and traveling wave branches are located at 5 ^ —L^. HoUinger et al. |5^, ^ have shown 
that r is a simple function of the amplitude and frequency of traveling waves and that the traveling wave 
branch can be divided into low-amplitude Soret and high-amplitude Rayleigh regimes. 

Many attempts have been made to reduce the governing-fluid dynamical equations to minimal mod- 
els which describe traveling waves. The most obvious approach is to extend the five-mode free-slip 
Veronis model ( [4.6D to include additional modes proportional to sin(A:z) |3^. However, the result- 
ing eight-mode model proves to be singular, because the Hopf bifurcation to traveling waves in binary 
fluid convection with free-slip boundary conditions is always degenerate |T9l pSl 27, 



|51| , as a con- 
sequence of (4.3). This means that the truncation (4.1), (|4.4|) is insufficient for even a qualitative de- 
scription of traveling waves in the full thermosolutal problem (|2.4| ). Numerous other models have ensued 
^ [sol , using other boundary conditions, additional field variables or different theo- 
retical approaches. Understanding the diverse aspects of traveling waves in binary fluid convection is an 
extremely challenging problem. 



6 Conclusions 

We have examined the well-known idealized thermosolutal problem from a variety of different perspectives. 

For infinite or large Prandtl number, the linear stability problem for temperature and concentration 
perturbations (T, C) is governed by a 2 x 2 matrix M whose entries depend linearly on the reduced 
Rayleigh number r, and whose eigenvalues provide the growth rates of perturbations to the motionless 
conductive state. 

We interpret the diagonal terms of this thermosolutal matrix as growth rates of two "pure" convec- 
tion problems driven exclusively by a thermal gradient or by a concentration gradient, which we term 
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the thermal and the solutal eigenvalues. Without coupling, the thermosolutal eigenvalues merely cross 
transversely as r is varied through the intersection point rjnt of the thermal and the solutal eigenvalues. 
Otherwise, the two eigenvalues of the coupled problem either undergo avoided crossing (the eigenvalues 
appear to deflect each other and remain real) or complex coalescence (the two eigenvalues join into a com- 
plex conjugate pair and then become real again). Which possibility is realized depends on the sign of the 
coupling: the product of off-diagonal terms. For the thermosolutal problem the coupling is proportional 
to, and has the same sign as, the separation parameter S. 

In the equivalent language of conic sections, the eigenvalues a + iuj and reduced Rayleigh number r 
satisfy a second-degree equation. Hence the sets (r, a) and (r, u) form hyperbolas, parabolas, or ellipses, 
according to the values of two invariants. The first is the discriminant of the matrix responsible for 
advection, which here is positive (except for the single value S = —1). The curves (r, cr) are thus 
hyperbolas and the curves {r,uj) ellipses. The second invariant A is, for the thermosolutal problem, 
proportional to the separation constant S. A = is the limiting case of a hyperbola consisting of two 
intersecting lines, or of an ellipse whose radii are zero. As A changes sign, the quadrants occupied by the 
hyperbola shift, and the ellipse becomes empty. 

Both of these equivalent formulations underline the organizing role played by the delimiting case 
S = 0. Most studies of binary fluid convection treat the codimension-two point ~ —L^, where 

the pitchfork and Hopf bifurcation curves meet at Tpp ^ = r^, = {1 — L)^^ , as a distinguished point in 
the (S'j r) plane, and expand around it. Our complementary point of view focuses on S" = 0, where the 
pure thermal and solutal eigenvalues intersect at Tint = 1 — L, as a different kind of distinguished point. 

Turning to the nonlinear problem, the minimal model of thermosolutal convection incorporating the 
lowest-order nonlinear effects was first derived by Veronis and has since been extensively studied. We 
find that the system of nonlinear equations satisfied by the steady states of the minimal model is of the 
special form: 

Mv^ E{vi,V2,-)v (6.1) 

where M is a matrix, v a vector, and E a scalar function of the components of v. The solution of such 
systems reduces to that of diagonalizing a matrix and solving a single nonlinear equation of one variable. 
The eigenvalues of M are the possible values of E. Its eigenvectors are used to reduce the number of 
arguments of E to one by rewriting all but one component of v as multiples of the remaining component. 

The five- variable system governing the steady states of the minimal thermosolutal model can be further 



reduced to a two-variable system in (T, C) of type (6.1) in which the scalar function E is proportional 



to the kinetic energy and the 2x2 matrix M bears a striking resemblance to the linear stability matrix 
M . This leads to a remarkable analogy between the linear stability problem and the nonlinear steady 
state problem. The energy E also undergoes avoided crossing or complex coalescence, again according 
to the sign of S. The curves (r, E) are hyperbolas. Complex coalescence for E must be interpreted as 
the disappearance of solution branches, rather than as the onset of oscillatory behavior. Quantitative 
results concerning the growth rates of perturbations to the conductive state can be translated to results 
concerning the kinetic energy of nonlinear steady states merely by transforming S —> LS and L ~> L^, 
where L is the Lewis number, the ratio of solutal to thermal diffusivities. 

For some results - the pitchfork bifurcation rpp - the analogy leaves the linear result unchanged. 
This reflects the bifurcation-theoretic fact that a change in sign in eigenvalue signals a bifurcation, i.e. 
an intersection between solution branches. For other results - the Hopf bifurcation Th - the analogy 
is valid, but not meaningful for the nonlinear problem. For yet other results - the point at which 
eigenvalues become complex r-|_ and the codimension-two point (S"*, r,^) - the nonlinear analogues are both 
significant and different from the linear phenomena. In particular, for negative S, the analogy indicates 
that there is a fundamental relationship between the onset of oscillation via the complex coalescence at 
r+ of two real eigenvalues in the linear problem, and the termination of two finite-amplitude solution 
branches via saddle-node bifurcation at f-|_ — for the nonlinear problem. As a corollary, there also 
exists a fundamental relationship between the merging of the pitchfork with the Hopf bifurcation in the 
codimension-two point (Bogdanov bifurcation) at (S'*,r*) for the linear problem, and the merging of 
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the pitchfork with the saddle-node bifurcation in a different kind of codiniension-two point (degenerate 
pitchfork) at {S^,f^) for the nonhnear problem. 

The relationship is a consequence of the exact analogy between the growth rates of the linear stability 
problem at infinite P and the energies of the nonlinear steady states of the minimal five-mode Veronis 
model both for the thermosolutal problem with imposed horizontally periodicity and free-slip vertical 
boundaries. We would not expect this analogy to persist when any of these assumptions are relaxed, 
i.e. when the boundary conditions are realistic, when the concentration gradient occurs via the Soret 
effect, when the amplitudes or the Prandtl number are moderate rather than infinitesimal or infinite. 
Yet, evidence |^, |2^, |2^, |5^, 57 suggests that the analogy must hold at least approximately, since 



the leading behavior of the two codimension-two points, S"* —L^ and S^, ~ —L^, and the domains of 
existence of the Hopf and saddle-node bifurcations continue to be related in approximately the same way 
as for the idealized thermosolutal problem. This indicates that the analogy between growth rates and 
energies of steady states could be a fundamental unifying feature of double-diffusive problems. 

Another provocative feature of binary fluid convection is the fact that the Rayleigh number is a single- 
valued function of any one of the following variables: growth rate, steady state energy, traveling wave 
amplitude and frequency, i.e. for any of these variables, each value is achieved at most once when varying 
r. We have explained this dependence for the growth rate and energy by the way in which r enters the 
problem via advection. HoUinger et al. [^6| |57| provide a related explanation invoking the reduction of 
the velocity field to one mode and the resulting simplification of the nonlinear terms. 

Finally, we have proposed a classification of eigenvalues and of steady states as primarily thermal or 
primarily solutal, based on their proximity to the eigenvalues or steady states of the pure thermal and pure 
solutal problems and on the relative proportions of solutal and thermal contributions to the buoyancy 
force. For many S values of interest, the coupling term is small in the vicinity of the bifurcations, and so 
the convection threshold is very close to that of pure thermal or pure solutal convection. The classification 
is particularly useful for the nonlinear problem for positive S. We have shown that the relatively abrupt 
transition between small amplitude and large amplitude convection called the Soret and Rayleigh regimes 
[ p2| |29| , ^ 1^ corresponds to the change in slope seen as a hyperbola adheres to first one 

and then the other asymptote as r is increased. This transition may be masked, because it occurs in a 
regime not corresponding to a real steady state, or muted, because it is too gradual. We have been able 
to give precise conditions under which the transition from Soret to Rayleigh regimes can be observed. 

We emphasize that our goal has not been to reproduce all of the spatio-temporal dynamics of double- 
diffusive convection, nor even to investigate the temporal dynamics of the five-mode Veronis model. 
Indeed, these goals have been admirably pursued and accomplished in previous research using other 
approaches, e.g. |l7[|^, |||^, |8[ H |3 • Nor can this purpose be accomplished by the idealized 



two-variable models we have investigated. Instead, our goal has been to extract certain universal large- 
scale features of double-diffusive convection in as simple a context as possible, and to re-examine these 
features in light of the avoided crossing/complex coalescence dichotemy and the linear/nonlinear analogy 
that we have put forth. 

Convection in binary fluids has previously provided a testbed for the discovery and realization of 
many fascinating phenomena in dynamical systems. Our hope is that these new perspectives continue 
this tradition. 
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Appendices 



A Finite Prandtl number model 



Although the results of section ^ of concerning the nonlinear steady states are independent of Prandtl 
number P, those of section || concerning the linear stability problem are derived by taking P to be infinite. 
In this Appendix, we describe our reduction of the thermosolutal linear stability problem to a 2 x 2 matrix 
in the case when P is finite. The key step in our interpretation is to decompose the velocity field into 
"thermal" and "solutal" velocity fields induced by the thermal and concentration gradients, with vertical 
components wt and wc, respectively. Referring to equations (3.1), the linearized equations governing 
this augmented set of fields are: 



dtq-^W^WT = Prk-'^dlf + Pq-'^V^WT 

dtC ^ WT + WC + Lq-^\/^T 
dtq-^V^wc = PSrk-^dlC + Pq-^V^wc 



Using the spatial and temporal dependence and notation defined in (p.2|)-(p.3|), (A.l) becomes: 
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(A.2) 



This 4x4 system has exactly the same eigenvectors and eigenvalues as system (3.4) with the additional 
eigenvector (T wt C wc) = (010 — 1) and eigenvalue fj, = —P. 

Note that the temperature and concentration fields are each advected by both the "thermal velocity" 
and the "solutal velocity", ft is this cross- ad vection which couples the thermal and solutal problem. 
Neglecting it leads to the decoupled thermal and solutal problems discussed below. 



A.l Thermal problem 



The upper left 2x2 submatrix of ( |A.2| ) describes the onset of thermal convection in a simple fluid of 
finite Prandtl number. This is perhaps the prototypical problem in hydrodynamic stability theory (e.g., 
[p^ Q 0)- The corresponding thermal eigenvalues satisfy 



T 
Wt 



-1 1 
Pr -P 



T 
Wt 



(A.3) 



P+1 



and are plotted in figure |l3[ 
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Pr 



(A.4) 



Figure shows that for r < — (P— 1)^/ (4P), eigenvalues are complex; perturbations to the conductive 
state oscillate a s the y decay. (In the terms introduced in the appendix, the discriminant of the r-dependent 
component of (A.3) is zero, so the curves of real values and imaginary parts of /LtT± are both parabolas.) 
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Figure 13: Eigenvalues of pure thermal problem for L — 0.01. Dotted and dashed curves shows [Lt for 
P = 1 and P = 10, respectively. Solid curve shows ctt, valid for P = oo, and also ht for P ~ 100, from 
which it is indistinguishable in this range. Coordinates of vertex are shown for P = 10. 



We consider only the regime in which the eigenvalues are real. Since we will only require the larger of 
the two eigenvalues, we will write /i-r = Mt+- For P large, ht becomes CTt of ( |3.8a ), as expected. The 
threshold of is 7't = 1, that of Ut. Indeed, as is well known, the threshold of convection is independent 
of P. The slope of fir at threshold is P/{P + 1), which also approaches that of ctt for P large. We will 
also require the normalized right and left eigenvectors corresponding to fiT- 
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(A.5b) 



where iV|, 



2^{{P - l)/2)2 + Pr. Thus 



(A.5c) 



where Mt is the pure thermal matrix in (A. 3). 



A. 2 Solutal problem 



The pure solutal problem, described by the lower right 2x2 submatrix of (A. 2) is completely analogous 
to the thermal case, with the inclusion of the Lewis number L and the separation parameter S: 
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The solutal eigenvalues, plotted in figure ^ are: 
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Figure 14: Eigenvalues of pure solutal problem for L — 0.01. and S — —0.1 (left) and 5 — 0.1 (right). 
Dotted and dashed curves shows [ic for P — \ and P — 10, respectively. Solid curve shows [ic for 
P = 100, which is indistinguishable from valid for P = oo. Coordinates of vertex are shown for 
P = 10. 



We again require only the larger of the two eigenvalues fic = f^c+, which approaches CTc for large P. 
The onset of convection occurs a.t rc = L/S, at which fic has slope PS/{P + L). We again consider only 
the regime in which the eigenvalues are real, i.e. ((L — P)/2)^ + PSr > 0, 



The corresponding normalized right and left solutal eigenvectors are: 



1 



1 



1 



where iV^ = 2y/{{P - L)/2y + PSr Thus 



C^McC" 



(A.8a) 
(A.8b) 

(A.8c) 



where Mc is the pure solutal matrix in (A. 6) 
A. 3 Thermosolutal coupling 



We now project the 4x4 thermosolutal problem onto the most unstable thermal and solutal mod es to 
form the 2x2 matrix which constitutes our approximation. We do so by multiplying the matrix of (|A.2| ) 
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by left and right eigenvectors as follows: 
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In ( [A.9b[ )-( |A!9c] ), Mt,Mc are the 2 x 2 pure thermal and solutal matrices and ,T^,C^ the corre- 
sponding left and right eigenvectors defined in (A. 5) and ([A.8|). B, T are the 2x2 oS'-diagonal submatrices 



in (A.9a). The calculation of the off-diagonal elements /3, 7 in ( |A.9c] ) is tedious but straightforward. Their 
product is: 

p2 2 

Pi = , (A.IO) 
4^((P - f)/2)2 + Pr)((P - L)/2)2 + PSr) 

This expression is not singular in the regime we consider here; the assumption that the eigenvalues 
jiT, ^J'C of the thermal and the solutal problems be real requires that both factors inside the square root 
be positive. The coupling ( A.1C| ), plotted in figure |l5| for P = 10, reduces to the far simpler coupling Sr'^ 
of (3.8c) in the limit of P = 00 and shares its salient feature: its sign is that of S, leading to avoided 
crossing if 5* > and complex coalescence if S* < 0. 



The eigenvalues of the 2x2 approximate matrix are: 




(A.ll) 



= MC+ 1 and f3i given by ( A^ ) , (|A.7| ) , and ( A.IO ) . In figures and |l^, we compare 



with ht = MT4 

results from: 

- the 2x2 approximate matrix ( |A.9c] ) whose eigenvalues /i± are given by ( A.ll ) for P = 10 (dashed 
curves) 

- the 3x3 exact matrix (3^) for P = 10 (dotted curves) 

- the 2x2 matrix ( |3.6|) whose eigenvalues a± are given by ( 3.11 ) for P — 00 (solid curves). 



Specifically, in figure [l^ we com pare the real parts of the eigenvalues ( A.ll| ) with those of ( 3.11 ) and 
with those of the two eigenvalues of (3.4) with largest real part. For the parameter values plotted, P = 10 
and S = ±0.1, all three expressions give very similar results. For both S = ±0.1, /x+ gives a slightly 
better approximation of the 3x3 eigenvalue than (t+, but cr_ is slightly more accurate than In figure 
p!?! , we compare the thresholds for pitchfork and Hopf bifurcations derived from these three matrices. The 
thresholds for the pitchfork bifurcations are so close as to be indistinguishable on the figure, whereas the 
Hopf bifurcation threshold is overestimated by the 2x2 approximate matrix. 

Aside from these quantitative comparisons, figures ^ and |l^ demonstrate that the exact finite Prandtl 
number linear stability problem and our reduced model both exhibit the essential qualitative features of 
thermosolutal convection: complex coalescence for S < and avoided crossing for S > Q; pitchfork, Hopf, 
and codimension-two bifurcations. 
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Figure 15: Product (3^ of off-diagonal terms. Dashed curves show expression ( A. 10 ) for P = 10, soUd 
curves show Sr'^ , vahd for P = oo. Positive values correspond to S" = 0.1, negative to S = —0.1. 
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Figure 16: Real parts of the eigenvalues /i± of the approximate 2x2 matrix (dashed curves) for P = 10 
and for S ~ —0.1 (left) and S — 0.1 (right). Shown for comparison are the real parts of the eigenvalues 
of the exact 3x3 matrix (dotted curves) for P = 10 and of the eigenvalues a± (solid curves) for P = oo. 
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Figure 17: Thresholds for pitchfork and Hopf bifurcations. Dashed curves are thresholds calculated 
from the approximate 2x2 matrix for P — 10. Dotted curves are thresholds calculated from the exact 
3x3 matrix for P = 10. Solid curves are thresholds for P = oo. The three thresholds for the pitchfork 
bifurcations cannot be distinguished; all lie on solid curves. 

B Conic sections and eigenvalues 

Consider a 2 x 2 matrix whose elements depend linearly on a parameter r: 

M = Mo+ rMi 



'a 13' 
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/3o ■ 




ai (3i 
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_ 7i Si 



(B.l) 



We wish to describe the dependence of the eigenvalues of M on r. The equation obeyed by the 
eigenvalues A(r) is: 

- TA + Det = (B.2) 

Here, and throughout this appendix, T refers to the trace and not to the temperature deviation; Det is 
the determinant of M. We have: 

T = {ai + 5i) r + ao + So 

= Tir + To (B.3) 



Det = (ai(5i - /3i7i) + (ao(5i + q;i(5o - /3o7i - /3i7o) ^ + "0(^0 - /3o7o 
= Deti + 2 Deti/2 r + Deto 



(B.4a) 
(B.4b) 



where Tq, Ti, Detp , Deti a re the traces and determinants of the matrices Mq and Mi in (BT) and Deti/2 
is defined by ( B.4a )-( B.4b ). The eigenvalues of M are complex where the discriminant Disc is negative, 
where 



Disc = T2-4Det = (T? - 4Deti) + 2(ToTi - 4Deti/2) r + (T^ - 4Deto) (B.5a) 
= Disci r'^ + 2Disci/2 r + Disco (B.5b) 
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Here Disco, Disci are the discriminants of Mq, Mi and Disci/2 is defined by (B.5a)-(B.5b). Whether and 
where Disc is negative is in turn seen from ( B.5b| ) to be determined by Disci and by 



A 



DisCj/2 ^ DisciDisco 
~ 4Disci 



(B.6) 



Disci and A are both invariants under translation and rotation of (r, A) of equation (B.2); see, e.g., 



Writing A = cr + iw, substituting (B.3) and (B.4b) into (B.2), and separating into real and imaginary 
parts, we obtain: 



Ticrr + Detir^ - Tqct + 2Deti/2r + Deto 



According to equation ( B.7b| ) 



either uj — 

or a-i(To + Tir) 



If a; = 0, then (B.7a) becomes: 

cr^ - Ticrr + Detir^ - Tocr + 2Deti/2r + Deto = 
If w 7^ 0, then substituting (|B.8b|) into (|B.7a| ) yields: 



n Disci o Disci/2 Disco 
uj^ H -r^ H —r H = 







(2cr - Tir - To) w = 



(B.Ta) 
(B.7b) 



(B.Sa) 
(B.8b) 

(B.9) 
(B.IO) 



The qualitative nature of the solutions to the second-degree equations (BJ) and( B.lO ) de pend s on 
the sign of Disci and A. If the discri minant Disci is positive, then the set (r, cr) satisfying (B^) is a 
hyperbola and the set (r, w) satisfying ( B.IO ) is an ellipse, and vice versa if Disci is negative. If Disci 
vanishes, then both sets are parabolas. In almost all the cases we shall study, Disci is positive and we 
shall assume this from now on. 



We rewrite equation ( |B.1C| ) as 
, Disci 

+ -r 
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where A is defined by (B^) and r,nid by ( B.lla )-( B.llb ). The sign of A is also crucial: if A > 0, then 
(B.ll) describes a degenerate ellipse containing no points: there are no complex eigenvalues. If A = 0, 
then the ellipse contains just the point (rniid,0). If A < 0, then ( B.ll ) describes an ellipse whose two 
semiaxes are r = r„ii(j of length \/— A and uj — of length -^Z— 4A/Disci; complex eigenvalues exist over 
the interval 

\r - n„id| < V-4A/Disci (B.12) 



At r — rmid, w attains its maximum value of \/— A and ( B.Sbj ) implies CT(rmid) = Cmid, where 
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(B.13) 



We now turn to the hyperbola described by (B.9) when Disci > 0; in particular we seek to characterize 
it by its asymptotes. The first three terms of (B.9) imply that the sum of the slopes of the asymptotes 
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is Ti and their product is Deti, i.e. the slopes are the eigenvalues Ai± = (Ti ± -\/Disci)/2 of Mi. These 
are real and distinct by the assumption Disci > 0. Some more algebra shows that (B.9) is equivalent to: 



(o- - CTinid - Ai+(r - r„iid)) (a- - (T„iid - Ai_(r - r^id)) = A 



(B.14) 



where rmid, A, and CTmid are defined in ( B.lla )-( B.llb ) and ( B.13| ) 



The asymptotes are the roots of the two factors in (B.14). They intersect at (rmid, cr,„id) and divide 



the plane into four quadrants. The magnitude of A measures the distance of closest approach of the 
two portions of the hyperbola; its sign determines which two of the four quadrants are occupied by the 
hyperbola. If A > 0, then one branch cr+ of the hyperbola lies above both asymptotes (in the sense of 
greater a) and the other branch tT_ lies below them. Each branch exists for all r. This is the situation 
called avoided crossing. When A = 0, the hyperbola is degenerate and consists precisely of the two 
intersecting asymptotic lines. In this case, the two branches f7+ and ct_ can be considered to either 
behave non-smoothly or to exchange identities at r = rmid- If A < 0, then both values a± lie between 
the asymptotes, i.e. they both lie above one asymptote and below the other. In this case, there are no 
real solutions to (B.14) in the range (B.12) surrounding rmid. At the endpoints of the interval in (B.12), 
the curves a± join, to be replaced within this interval by the single linear segment (B.St). This is the 
situation we call complex coalescence. 



The phenomenon of avoided crossing can be quantified by differentiating (B.14) implicitly with respect 



to r: 



(cr' - Ai+)((T - (Tmid - Ai_(r - Tmid)) + (ct' - Ai_ ) (cr - (Jinid - Ai+ (r - r^id) ) = (B.15a) 



a"(2((T-(Tmid)-(Ai_+Ai+)(r-rmid)) + 2((t' - Ai_)((t' - Ai+) = 
(t"'(2((7- CTmid) -(Ai_+Ai+)(r- rmid)) + 3(t"(2(7' - (Ai_ + Ai+)) = 



(B.15b) 
(B.15c) 



and then evaluating cr, a', a", and a'" successively at rmid: 
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(B.16a) 
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(B.16c) 
(B.16d) 



We see from ( B.16 ) that the change in slope undergone by cr± at rmid corresponds to an extremum in 
cr^, whose magnitude measures the abruptness of the change. 

Another fact which we shall use is that a line in the (r, cr) plane which is parallel, but not equal, to 
one of the asymptotes intersects the hyperbola in exactly one point, 



We briefly discuss the exceptional case Disci = 0. Equation (B.IC) for the imaginary part of the 
eigenvalues becomes: 
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(B.17) 



Equation (B.9) for the real part of the eigenvalues becomes: 
= (cr-^r-^l + f 2Deti , 



Ti To 

= 10" ^r — 

2 2 
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TqTi 
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(B.18) 



Both (B.17) and ( B.18| ) describe parabolas. The parabola of ( B.1'1) is o riented along the the r-axis, while 
the axis of ( B.18 ) is the line a — (Tir + To)/2. If Ti — 0, then ( fB.T^ ) is also oriented along the r-axis. 
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The two parabolas are oriented in opposite directions: ( |B.17 ) opens towards positive values of r if Disci/2 
is negative and vice versa for (B.18). The vertex of both parabolas is located at r = — Disco/(2Disci/2)i 
with w = for ( |B.17D and a = (-DiscoTi/(2Disci/2) + Tq)/2 for ( |B.18D . 



Our treatment of binary fluid convection leads to two matrices of type (B.l), one whose eigenvalues 
a govern the linear stability of the conductive state and the other whose eigenvalues E are the kinetic 
energy of nonlinear steady states. L is the Lewis number, which is necessarily positive and usually small, 
iS" is the separation parameter, which may have either sign, and r is the reduced Rayleigh number. In 
tables 1 and 2, we give the quantities we have defined above for each of these two matrices. 



Growth rate matrix 
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